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ITERATIVE  ABEL  INVERSION  OF  OPTICALLY  THICK 
CYLINDRICALLY  SYMMETRIC  RADIATION  SOURCES 


The  transverse  monochromatic  radiance  and  transmittance 
profiles  for  a  cylindrically  symmetric  radiation  source  are 
related  to  the  radial  emission  and  absorption  coefficients  by 
the  radiative  transfer  equations ^ 


R 


S(z)  =  2  I  J(r)  cosh  G(z,  r) 


r  dr 


(r2-s2,1/2 


(1) 


and 


R 


-lnT(z)  =  2  K  (r) 


/ 


r  dr 


,  2  2.1/2 

(r  -  z  ) 


■  =  2  G(z,  R) 


(2) 


where  S(z)  =  N(z)/t^(z);  N(z),  t(z),  and  z  are,  respectively, 
the  transverse  radiance,  transmittance,  and  coordinate;  J(r), 
K(r),  and  r  are,  respectively,  the  radial  emission  coefficient, 
absorption  coefficient,  and  coordinate;  R  is  the  source  radius; 


and 


These  equations  are  valid  for  arbitrary  optical  thickness  of  the 
source.  The  problem  addressed  is  the  solution  for  J(r)  and  K(r) 
from  N(z)  and  t(z)  by  inversion  of  Eqs.  (1)  and  (2).  No  considera¬ 
tion  is  given  here  to  the  problem  of  random  or  bias  error  propa¬ 
gation  in  the  inversion. 

Regardless  of  optical  depth,  the  transmittance  equation  is  an 
Abel  integral,  and  K(r)  can  be  found  immediately  from 
K(r)  =«Vt-ln  t(z)]  where  is  the  Abel  transformation  operator 


.V[f(z)] 


(4) 


For  optically  thin  sources,  that  is, for  RK(r)  «  1  for  all  r, 
cosh  G(z,  r)  **■  1  and  the  radiance  equation  also  reduces  to  an 
Abel  integral  equation.  Then,  Eq.  (4)  can  be  used  to  obtain 
J(r)  =  .V  [S(z) ]  .  If  the  source  is  not  optically  thin,  an  iterative 
procedure  can  be  used  to  obtain  J(r). 

Elder  et  al.  ^  proposed  the  following  iterative  solution 


J(r)  =  JQ(r)  +  J2(r)  +  J2(r)  +  ... 


where 


J0(r)  =  ,V[S(z)] 


(5) 


and 


J„00  =  .v 


XN. 


(r)  Lcosh  G(z,  r)-l  J 


r  dr 

(r2-*2)1'2 
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The  summation  is  truncated  when  Jn(r)  is  less  then  some 
preassigned  minimum.  The  authors  claimed  rapid  convergence 
with  this  method  and  quoted  a  4%  accuracy  with  only  two  iterations 
for  values  of  the  centerline  transmittance  t(z  =  0)  as  small  as 

~0. 2.  Unfortunately,  although  the  convergence  is  rapid,  it  is 
not  to  the  correct  solution.  The  flaw  in  the  scheme  is  the 
implicit  correction  of  J(r)  based  on  an  Abel  inversion  of  the 
difference  between  the  previous  iteration  result  for  S(z)  and 
the  effective  thin  source  value  of  S(z),  that  is,  the  value  of  S(z) 
computed  from  Eq.  (l)with  J(r)  taken  as  Jg(r).  Intuitively, 
the  correction  should  be  based  on  an  Abel  inversion  of  the 
difference  between  the  previous  iteration  for  S(z)  and  the 
known  profile  S(z).  The  iteration  scheme  based  on  this  reason¬ 
ing  is 


Jn<r>  =  Jn-l<r> 


XV 

Jj- 


(r)  cosh  G(z ,  r) 


r  dr 


.  2  21/2 
(r  -z  ) 


-S(z) 


where 


(6) 


J0(r)  =^[S(z)]. 

d 

Here,  Jn(r)  is  the  nth  iteration  result  for  J(r)  and  not  a  correction 
term  as  in  Eq.  (5). 
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Inversions  with  Eqs.  (5)  and  (6)  were  carried  out  for 
simple  example  test  cases  in  which  both  K  and  J  are 
constants.  From  Eqs.  (1) and  (2),  the  transverse  input  profiles 
are  then  found  to  be 

S(z)  =  sinh  [K  (R2  -  z2)  1/2 ] 

and  (7) 

-lnT(z)  =  2K(R2-z2)1/2. 

R  and  J  were  set  to  unity,  and  K  was  selected  to  give  the  center 
line  transmittance  t(z  =  0)  =  0.9,  0.8,...,  0.1.  Inversions  of  these 
profiles  were  then  carried  out  with  the  two  methods.  Iteration  was 
continued  until  J(r  =  0)  differed  by  less  than  1%  between 
iterations.  A  numerical  quadrature  formula  based  on  the  method 
of  Barr^'  ^  was  used  for  both  the  iterative  calculation  of  S(z) 
and  the  Abel  transformation.  The  final  J  retrieved  from  the 
two  inversion  algorithms  are  shown  in  Fig.  1.  It  is  evident 
that  the  method  of  Elder  et  al.  does  not  converge  to  the  correct 
solution,  whereas  the  present  method  does.  The  rate  of  convergence 
is  the  same  for  the  two  methods  as  is  indicated  at  the  top  of  the 
figure  by  the  number  of  iterations  required  to  achieve  the  conver¬ 
gence  criterion. 
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ELDER  et  al.  ALGORITHM,  EQ  151 


i 


Extensive  verification  and  use  of  the  present  inversion 
algorithm  has  been  made  within  the  context  of  band  model 

(3) 

radiation  formulations  for  axisymmetric  rocket  plume  sources. '  ’ 
In  that  application,  the  only  essential  differences  from  the  present 
monochromatic  application  are  that  the  kernel  function  cosh  G(z,  r) 
is  replaced  by  a  different  (and  more  complicated)  function,  the 
radiance  equation  is  formulated  for  N(z)  rather  than  S(z),  and 
the  absorptance  equation  is  no  longer  an  Abel  integral  equation 
regardless  of  optical  thickness.  The  significant  consequence  of 
the  latter  difference  is  that  in  the  band  model  application,  simul¬ 
taneous  iterative  solutions  must  be  carried  out  for  both  J(r)  and 
K(r) .  The  method  was  shown  to  be  stable  and  accurate  for  a  wide 
range  of  values  and  profiles  for  J(r)  and  K(r). 
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pereonnel  in  dealing  with  the  many  problema  encountered  in  the  nation'a  rapidly 
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Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia¬ 
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